The persistence of viscous effects in the overlap region, 
and the mean velocity in turbulent pipe and 
channel flows 

Katepalli R. Sreenivasan and Anupam Sahay 

Mason Laboratory, Yale University 
New Haven, CT 06520-8286 



Abstract 



We expand on our previous argument ( fjreenivasan 1987 ) that important elements of the dynamics of 



wall-bounded flows reside at the wall-normal position y p corresponding to the peak of the Reynolds 
shear stress. Specializing to pipe and channel flows, we show that the mean momentum balance in 
the neighborhood of y p is distinct in character from those in the classical inner and outer layers. We 
revisit empirical data to confirm that y p = 0((hu /U*) 1 / 2 ) and show that, in a neighborhood of order 

1/2 

Rj around y p , only the viscous effects balance pressure-gradient terms. Here, R„ = hU*/v, h is 
the pipe radius or channel half-width, v is the kinematic viscosity of the fluid and U* is the friction 
velocity. This observation provides a mechanism by which viscous effects play an important role in 
regions traditionally thought to be inviscid or inertial; in particular, it throws doubt on the validity of 
the classical matching principle. Even so, it is shown that the classical semi-logarithmic behavior for 
the mean velocity distribution can be a valid approximation. It is argued that the recently advanced 
power-law profiles possess a rich underlying structure, and could be good approximations to the data 
over an extended region (but they too are unlikely to be exact). 

1 Introduction 

Essentially all the important notions in wall-bounded flows are cast in terms of two length 
scales: the inner or viscous scale v/U*, where v is the kinematic viscosity of the fluid 
and £/* is the friction velocity (= y/r w /p, t w and p being the wall shear stress and fluid 
density, respectively), and the outer scale h, where h is the pipe radius, channel half-width 
or boundary layer thickness. The ratio of the two scales is the Reynolds number i?* = hU*/v. 
Since the appropriate asymptotics correspond to the limit R* — > oo, one may expect that 
the problem has elements of singular perturbation. We shall use the standard notation that 
U + = U /£/* and y + = yU^/v, y being the normal distance from the wall. It is traditionally 
thought that the viscous effects are important up to a y + of about 30 and that, within this 
region, the outer length scale h is unimportant. In the bulk of the flow excluding this viscous 
region, it is thought that v is unimportant and the characteristic length is h. This view has 
been quite successful in organizing various experimental data (see, e.g., Sreenivasan 
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Dussauge et al. 1996 ), though it has been recognized for some time (e.g., Rao et al. 197l| ) 



that the interaction between the two scales is the key to the flow structure. 



In an earlier paper ( Sreenivasan 1987 ), we noted that a proper understanding of the 
boundary layer structure requires greater emphasis on wall-normal position where the Reynolds 
shear stress peaks. This peak position scales as the geometric mean of the inner and outer 



length scales. It is well known (e.g., Drazin fc Reid 1981) that, in linear and early stages of 



nonlinear instability in boundary layers and channel flows, the position of the peak Reynolds 
shear stress coincides with the critical layer. This observation was inverted in |Sreenivasan 



(1987 ) to suggest that the position of the peak Reynolds shear stress, y p , in turbulent wall 
flows plays something of the same role as that of the critical layer in unstable wall flows. It 
was pointed out that, just as for critical layer in unstable wall flows, the mean velocity at 
y p is approximately a constant fraction of the freestream or centerline velocity, the fraction 
being about 0.65. Other analogies between the critical layer in the unstable state and the 
'critical layer' in the turbulent state were also cited. 

The critical layer in unstable boundary layers is the seat of perturbation vorticity which 
undergoes amplification when the Reynolds number exceeds a certain threshold. The next 
stages of the perturbation development involve the onset of three dimensionalities and, 
eventually, of turbulence itself. If the turbulent vorticity in the boundary layer can be 
caricaturized as a vortex sheet, it was thought by analogy that its seat would be y p . By 
interaction with image vorticity (invoked to mimic the presence of the wall) , this hypothetical 
vortex sheet located at y p gets lumped first into two-dimensional rolls and, eventually, into 
horse-shoe shaped vortices. These latter structures are comparable in several respects to 
those found from visualization studies of the boundary layer (e.g., Head & Bandyopadhyay 



1981 ). It was further argued, albeit with less certainty, that the same picture can explain 
aspects of the structure in the wall region, for example the (noisy) spanwise periodicity of 
streaks ( |Kline et al. 1967| ). For quantitative details and comparisons with data, one should 
consult preenivasan (1987| ) . 

Whatever the detailed objections to the physical content of the model and however 
preliminary the attempt, it appeared that a self-contained picture of the boundary layer 
could be developed on that basis. Unlike in the predecessor paper ( preenivasan 1987 ) which 



focused on the large-scale structural elements of the boundary layer, we shall examine here 
the mean velocity distribution to reiterate, in quite a different way, the importance of the 
'critical layer'. This seems to be an especially timely goal because of the renewed interest and 



recent controversy surrounding the mean velocity distribution in pipe flows (e.g., Barenblatt 
fc Chorin 1996j |Zagarola fc Smits 1991 ). 



In the classical picture, one assumes the existence of a common region of validity of the 
outer and inner solutions of dU + {R 1f ,y + ) / dy + in the limit i?* — > oo,y + — > oo, and puts 
forth asymptotic arguments to obtain 

U+(R*,y+) = (l/K)lny+ + B, (1) 

where k and B are empirical constants presumed to be independent of R*. This is the 



celebrated log-law, which occupies a central place in the turbulence literature (e.g., Coles fe 
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Hirst 1969, Monin &; Yaglom 1971 , Tennekes fc Lumley 1972] ). It is generally thought that 
y + ~ 30 and y + ~ are the lower and upper limits of the logarithmic profile and the 

Karman constant k ~ 0.4 and B ~ 5.5 (see, for example, |Coles fc Hirst 1969 ). There is 
still some uncertainty about these constants: recent high-Reynolds-number measurements 
m pipe flows ( |Zagarola 19~96| ) yield k = 0.44 and B = 6.1. 

There exist alternative formulations for the overlap region ( [Long fc Chen 1981| , [Baren 



blatt 1993| , |Barenblatt fc Chorin 1996; , |George et al. 1996; , |George fc Castillo 1996| ). Here, 
we shall restrict attention to Barenblatt's formulation. Its primary contention is that the 
limit of small-viscosity (or high Reynolds number) is singular — as is common in second- 
order phase transitions in condensed matter ( Domb &: Green 1976] ) and also, perhaps, in 
Kolmogorov turbulence (Monin &; Yaglom 1975) — and so the viscous effects never disappear 
in the overlap region. This imperils the classical matching argument and the orthodox view 
that the log-law is exact in the infinite Reynolds number limit. Specifically, note that di- 
mensional considerations allow us to write the velocity distribution in an intermediate layer 
in the form 

y + {dU + /dy + ) = 4>(y + ,R*), (2) 

where if) is an unknown function of its arguments. In the classical picture, ip is thought to 
asymptote to a constant, say 1/ac, as the arguments of -0, namely i?* and y + , assume large 
values. Integration then yields the log-law. On the other hand, suppose that 



y + {dU + /dy + ) = {l/K){y + Y 



(3) 



a being some positive constant. This leads to a power-law for the mean velocity distribution. 
In particular, Barenblatt (1993) and Barenblatt & Chorin (1996) predict the asymptotic 
nature of U + in two regions of the flow — the classical overlap region (say Bl) and a region 
further out towards the center (say B2). In the power-law paradigm, Bl + B2 together form 
the overlap region. The specific predictions are the following: (a) In Bl, U + is tangent to 
the classical logarithmic profile to which it remains close but from which it always remains 
distinct, (b) In B2, the power law can be approximated properly by a logarithmic function 
similar to the classical log-law but with a slope that is approximately ^fe times that of the 
classical value. It may be thought that this latter prediction does not contradict the classical 
log-law because the domains of the two logarithmic regions are disjoint. However, because 
B2 would be a part of the outer region in the classical picture, one may consider that a 
conflict does exist here as well. 

The purpose of this paper is to examine the nature of the mean velocity distribution 
briefly, emphasizing along the way two significant qualitative issues: (a) The viscous effects 
are important in a region of pipe and channel flows that is traditionally thought to be 
inviscid and that, in fact, the balance there is between viscous and pressure gradient effects. 



Following Long fc Chen (1981 ), we might call this 'critical' region a mesolayer (although we 
do not necessarily subscribe to all the implications of that work). The existence of such a 
mesolayer gives a new twist to the dynamics of the boundary layer, but the degree to which 
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the classical picture needs modification is not yet clear (see section 4). (b) The importance 
of viscosity in the mesolayer offers a key to the regeneration mechanism of the boundary 
layer. The discussion here will be necessarily brief, and more details will be published 
elsewhere (Sahay & Sreenivasan T996| ). 



2 The background 

2.1 The wall-normal position of the peak of the Reynolds shear stress 

Central to the present arguments is the manner in which the Reynolds or turbulent shear 
stress, t = —{uv), is distributed in the flow; here u and v are velocity fluctuations in the 
streamwise coordinate x and the wall-normal coordinate y, respectively, and (•) denotes a 
suitable average. Its behavior in pipe and channel flows is shown qualitatively in Fig. 1. The 
quantity r + = t/U% increases from its zero value at the wall, apparently like y 3 for small 
y, and rapidly reaches about half the maximum value at a y + of about 12. It continues to 
increase further to reach a maximum value at y p . The maximum value equals t w in the limit 
-R* — > oo, but falls short of t w at any finite R*, by an amount, say, A(i?*). The Reynolds 
stress r decreases beyond y p , and reaches zero in the freestream of the boundary layer, and 
on the center line of the pipe or channel. 
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Figure 1: A schematic of the turbulent shear stress t + profile in channel and pipe flows, showing 
three distinct regions. The buffer layer is interposed between regions I and II. 
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Of special importance is the position y p . Its leading or der variation has been obtained 
empirically by Long fc Chen (1981 ) and Sreenivasan (1987 ), who have shown that 



Vt = A^ 1/2 , (4) 

where A = 1.87 and 2, respectively. Although the two prefactors are somewhat different, 
they agree on the principal result that the peak of the Reynolds shear stress occurs at a y + 
that increases as R l J 2 . We revisit this issue here. Figure 2a shows plots of r + vs y + /R^J 2 

1 /2 ! 

for a range of Reynolds numbers. The Rj variation of y p appears to be a good leading 
order approximation. There is some correction to this scaling at the lowest R*, which we 
shall examine subsequently. For the present, we have ignored the low i?* data in estimating 
A. On the basis of Fig. 2a, we take R* < 500 as an operational definition of the low Reynolds 
number. Figure 2b is an expanded plot near the peak of r + . An accurate determination 
of y+ for large R* is difficult because the peak is rather flat (and becomes more so with 
increasing R*) and because there is much scatter in the data. Mindful of these uncertainties 
we estimate that A = 1.8 ± 0.2. 

! 1 /2 

The point to emphasize is that, for all but the very low R*, y^ = 0(Rj ) lies well 
within the classical logarithmic region (30 < y + < 0.15-R*). We shall now discuss its role in 
determining the distribution of the mean velocity in boundary layer flows. The discussion is 
specialized for analytical convenience to plane channel flow and axisymmetric pipe flow. The 
simplicity to be gained is that all the terms in the mean momentum equation are independent 
of the streamwise direction. 



2.2 The basic physical idea 

The exact mean momentum equation is given by 

_dP + dr +u ^V_ _ Q ~ 
dx dy dy 2 

where P is the mean pressure and U is the mean velocity depending only on y. The Reynolds 
shear stress term appears in the equation as an unknown. We now make the obvious point 
that, at the position at which the turbulent shear stress r is a maximum, i.e., at y p , the 
pressure gradient terms are balanced only by viscous terms; the Reynolds stress terms are 
entirely absent because the quantity that appears in the momentum equation is the Reynolds 
shear stress gradient, not the Reynolds shear stress itself. We have already seen that y p 
resides in the part of the boundary layer traditionally thought to be independent of viscosity, 
or purely inertial. This means that, within a region in pipes and channels that has been 
thought to be inertial, there exists a neighborhood within which only viscous terms are capable 
of balancing the pressure gradient terms nearly entirely. 

How large is this neighborhood? To answer this question, we plot in Fig. 3 the ratio 
of the Reynolds shear stress gradient term to the viscous term in the vicinity of y^. Using 
experimental data for the purpose would have generated much scatter, so our estimates are 



5 



T 



+ 



0.8 



0.6 



0.4 



0.2 



o 



o ° o_Ooo* 




100 



-+ 



0.95 




0.85 







522 


— i — 


715 


— -x- — 


851 


—-m — 


916 


a 


1020 


— 


1177 


— o — 


1430 




2275 


A- 


4324 


A 


. 7309 


— v 


£8530 


-------- 



Figure 2: Plots of the turbulent shear stress t + as a function of y + / 'R\' 2 (a) across the channel and 
(b) near its peak for high Reynolds number experiments. The sources for the experimental data are 
Antonia et al. 1992 (channel, R* = 256, 916), Comte-Bellot 1963 (channel, i?„ = 4324, 7309), Kim 
et al. 1987 (channel DNS, i?* = 180), Laufer 1950 (channel, = 522, 1177, 2275), Laufer 1954 
(pipe, R* = 8530), Sirovich et al. 1991 (channel DNS, R* = 125), Wei & Willmarth 1989 (channel, 
i?* = 715, 1020), and Zagarola 1996 (pipe, i?* = 851, 1430). The shear stress has been obtained by 
the numerical differentiation of the measured velocity profile using Eq. (^|) of section 3. Zagarola's 
data for higher Reynolds numbers could not be used because the mean velocity data have not been 
measured close enough to the wall. 
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Figure 3: The ratio of the Reynolds shear stress gradient to the viscous stress gradient in region 
II around y p . The four curves are obtained for different fits to the mean velocity data, as explained 
at the bottom of the figure. The ratio for the Barenblatt profile is given by (1 + £) Q — 1 where 
a = 3/(21ni?e). 



based on various fits to the data: the classical logarithmic profile and various power-law 
profiles recommended by [Barenblatt (1993| ), These estimates suggest that the turbulent 
stress gradient term is of the order of a tenth of the viscous stress gradient term within the 
region —0.1 < £ < 0.1, where £ = (y + — y+)/y+ . In an order of magnitude sense, the region 

1 /2 

around y+ where viscous terms overwhelm turbulence terms has itself a width 0{Rj ). 

3 The momentum equation 

Integrating Eq. (||) and applying the boundary condition at y = to obtain the constant of 
integration, and that at h to eliminate the pressure gradient term, one obtains 

dU+ ^ + 1 y + ™ 
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It appears natural that we should use Eq. @ to understand the mean velocity distribution 
by approximating r + in the form of a double expansion around y^ and about the supremum 
value of unity (attained in the infinite Reynolds number limit). We have accordingly per- 
formed a local analysis (local in y + ), in the limit 12* — > oo, of Eq. (|fj) in regions I and II and 
III shown in the typical r + profile of Fig. 1. The regions are defined as domains of validity 
of the asymptotic expansion of dU + /dy + to the leading order in the following limits: 

12* —> oo, y + = 0(1) (region I), 

12* -> oo, y + = 0(y+) (region II), 

12* —* oo, y + = 0(12*) (region III). 

I is the classical viscous region, and III can be thought to be the classical outer region. 
Viscous terms are significant in both I and II, but the classical buffer region, in which these 
terms are small in relative magnitude, is interposed between the two regions. This makes 
the two regions distinct. Here we shall focus on II nearly entirely. The analysis of regions I 
and III in the same spirit yields some significantly new results ( [5ahay Sreenivasan 
which will be summarized only as needed. 

3.1 Analysis of region II 

We can write the Taylor series expansion of r + about y+ as 

oo 

r + {y + )=r+ ax + Y J fn(R*)C 

n=2 

where f = (y+ - y+)/y£, n\f n /y+ n = [d n r+/d n y + ](y+=y+), /„ = o(Rj) where 7 is any 
positive number, and T+ ax is the maximum value of r + for any given 12*. Clearly, r^ ax < 1, 
the equality holding only at infinite Reynolds number. For any finite Reynolds number, we 
write T+ aa , = 1 — A (12*) (see Figs. 1 and 2), where A = o(l). We may thus rewrite the 
above equation as 

00 

r+(y+) = [1 - A(12*)] + ]T / n (JJ.)f\ (7) 



n=2 



Putting (0) into @ we get 



rUT+ V + V + °° 



dy+ R* 12* 

n=2 



If we take dU + /dy + to be positive and monotonic with respect to y + (which certainly 
seems to be case empirically, although a theoretical proof is lacking), we may argue that 
A = 0(yp / R*): positivity implies A12*/y+ = o(RJ) for any positive 7 and monotonicity 
implies y+/A12* = o(Rj). 
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The result that A = 0(yp /R*), when used in conjunction with the estimate of y+ (Eq. 

(4)), yields A = 0{R* ), say ai?* , where a is a constant. (This latter relation is verified 
independently also by experiment.) Using these estimates of the leading order of and A 
we get, from Eq. <M), 



y ±t* it* it*. ii—o 



dll + a A A '^-l/2s 

L x t <,j -r 

n=2 

Making use of the relation y + = yt{l + £) = XR\ (1 +() we can rearrange (||) to obtain 



_ (o_A)A + (o-2A)A + 1/2 £ 
dy + y+ y + r ^ 



3.2 The logarithmic 'law' 

An exact logarithmic profile would obtain if, in a nonzero neighborhood of £ = 0, the first 
term in Eq. (|io| ) is dominant in the limit R* — > cc,y + = 0(y^"). The second term is of 
the order of the first term and thus can be neglected only if it is identically zero, which 
requires the exact equality of a and 2A. This assessment requires a numerical estimate of 
the constant a, which we examine in Fig. 4. 

Here again it is difficult to determine the coefficient accurately owing to finite- Reynolds- 
number corrections. Ideally, one would consider data in intervals like [R*, oo) for increasing 
R* and look for convergence of the least-square estimates of the coefficient a. The inset in 
Fig. 4 shows our attempt at implementing this algorithm within the constraints of finite 
maximum R* and small number of data points. Least-square estimates of a are obtained 
using n data points at the largest available Reynolds numbers. We have taken into account 
only the data for R* > 500. We do not consider smaller than five points (i.e., we stop at 
n = 5) since the fit for a smaller value of n will lead to unacceptable statistical uncertainties. 
The resulting estimate for a is 3.1 ± 0.1 (see inset). 

The uncertainty in the numerical estimates of a and A allow for the possibility that 
a — 2A = 0, although the use of their mean values makes the second term nonzero. It is clear 
that empirical estimates will not settle the issue of exact equality of the two numbers. An 
added ignorance factor is the lack of knowledge of the asymptotic behavior of the functions 
fn{R*)- Fo r the logarithmic law to exist, the f n (derivatives of r + at y p ) must be o(A). 

In spite of the ambiguity regarding the existence of an exact logarithm law in region II, it 
is true that there exists a nonzero neighborhood around y^ in which a logarithmic variation 
is a good approximation to the true velocity profile. In the extended limit R* — > oo,y + = 
0(yt),£ — ► 0, the leading term representation of ( |IoD is 

dU+ _ (a-A)A 
dy+ y+ 



9 




0.02 0.04 0.06 0.08 0.1 

Figure 4: Reynolds number variation of the leading order correction to the peak value of r + . 
The convergence of the estimate of a was checked by varying the number of data-points (n points 
corresponding to the largest i?*) used for the least-square fit. The variation of a with n is shown in 
the inset plot. The error bars represent 68% confidence intervals. 



just as required for a logarithmic behavior. The log-law constant 1/k will then be 2.34 ±0.8. 
Despite the large uncertainty, the mean is very close to the traditionally accepted value. A 
detailed analysis of region III of Fig. | flgahay fc Sreenivasan 199(f ) shows, however, that the 
outer edge of this logarithmic region is o(l). 



3.3 Local structure of the power-law profile 



The expression for the profile is ( Barenblatt 1993 , Barenblatt fc Chorin 1996 ) 

3 



u- 



a 



2 In Re ' 



(12) 



where Re is the Reynolds number based on the pipe diameter and the average velocity 
across the pipe cross-section. The factor v3 in (3 is an aesthetic choice of the originator of 
the equation. The basis for the choice of In Re in a is that In Re is insensitive to the precise 
definition of Re (see, e.g., Barenblatt fc Goldenfeld 1995|). The rel ationship between Re and 
R* is implicit in ( |l2|) and has been shown to be ( |Barcnblatt 1993| ) 



R* 



e 3 / 2a 2 a a(l + a)(2 + a) 
\/3 + 5a 



(13) 
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Consider the limit process /?* - 
Ua)/yi ■ We can write Eq. (|l2|) as 



0(lla) where y+ is arbitrary. Let £ = (y H 



U + = /3e 



olny+ aln(l+C) 



(14) 



Since C is 0(1) and a is 0((lni2e) x ) the second exponential can be expanded to yield, after 
some rearrangement, 

(aPyf) lny + + (3y+ a (l - lny+ a ) + 0((a ln(l + C)) 2 )- (15) 



Thus in a small neighborhood around y+ the mean velocity is like a logarithm. The asymp- 
totic forms of the slope (A = a{3y£ a ) and the constant (C = /3y+ Q (l — lny+ Q )) will depend 
upon the order and magnitude of y+. If we take y+ = aRj then it follows that 



A 



C 



,3 7 /2 



Vs Q( ln2R * 

2 AnR* 



lni?* + 0(ln 2 i?*) 



(16) 
(17) 



where In 2 R* = ln(lni?^). Putting 7=1 the logarithmic variation predicted for region III 
( [Barenblatt fc Chorin 1996 ) is obtained. The slope of this logarithmic law is indeed y/e 
times larger than the universal log-law. 

It is of interest to examine the local structure of Eq. ( |T2] ) near y p . The local logarithmic 
approximation (|i~5| ) will hold with y+ replaced by y p . To calculate the asymptotic forms 
of the slope and the constant (labeled A p and C p respectively) we need an expression for 



Vp which is consistent with the power-law expression for U + , Eq. (12). This can be 



done by substituting ( |12D in the mean momentum equation 
which dT + /dy + = 0. These manipulations yield 



(1 — a)a[3R* 



l/(2-a) 



and solving for y + = y+ at 



(18) 



In the limit iL 



oo, Eq. ( |l8| ) can be written as 

1/2 



Up 



R* 



1/23/8 




n=l rra=0 



(In 2 #,) m 
(lni?*)" 



(19) 



where P mn are constants. The asymptotic forms of A p and C p can easily be obtained. 

The form of U + at y p corresponding to the power-law ( |i~2|) can be derived in a rational 
way from the mean momentum equation, along the lines of the analysis in section 3.1. We 
take into account (possible) corrections to the leading order of y p and assume 
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1/2 



n=l m=0 



(In 2 R,) rn 
(lnRA n 



(20) 
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Figure 5: The difference between A and the ieading order behavior 3.1i?» ' , plotted against 
In 2 i?*/(i?y 2 In i?*), in order to determine the next-order term in Eq. (|2l|). Using the best fit to the 
data (except for those at the lowest three Reynolds numbers) yields the relation A = 3.1i?* (1 + 
0.3 In 2 i?*/lni?„). The symbols key is identical to that of Fig. 4. 



in accordance with (|l9|). It is prudent to resolve A to the same degree as y + , hence we take 

(In 2 R*) r 



A = aR* 



-1/2 



oo n 

1 + W a 



7 „ 7 > u " n!ir 

n=l m=0 



(lnii,)" 
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The logarithmic 'correction terms' to the leading order dependence of A on R* are not 
unrealistic. For example, returning to Fig. [|, it is seen that the R* dependence is good 
only at very large i?*. In Fig. ||, we show the difference A — aR* 1//2 , which itself behaves to 
the leading order as In 2 R*/{R^ 2 lni?*), just as supposed in the expansion (pTl) . 



(In 2 IQ 



Using Eqs. (20), (|2l|) and the relation y + = y + (l + (), we obtain an expression similar 

oo n 

1 + E E <4 



to (10) for the velocity derivate 
dU+ 



dy^ 



aX , „. 

— i + C 

y + 



A 2 (i + C)' 
y + 



n=l m=0 

OO 'I 



* + EE^ 



{lnR,) n 
(In 2 R, 



where a' 



n=l m=0 

and A' are constants related to a m » 



(InR*) 



and X r , 



+ E/nC" (22) 
n=2 

In the extended limit R* — » 



oo, y + = 0(yp), C — » the leading term of ( ]22| ) is of the same form as that of d£/ + /dy + for 
the Barenblatt profile. 
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In summary, it appears that the Barenblatt-Chorin profile is consistent with Eq. @ 
when the double series expansions © and (H) are used for y+ and A. 

4 Conclusions 

The primary qualitative point of this work is that there exists a mechanism for viscous 
effects to spread to the classical overlap region. The existence of such a mechanism will in 
principle prevent classical matching, so that the leading order 'inner' and 'outer' expansions 
are technically insufficient to construct a uniformly valid approximation of the mean velocity 
profile. 

What, then, can be said about the classical log-law? First, it must be noted that the 
velocity change across the 'mesolayer' is not of the order unity, unlike in most boundary 
layer problems, which may make its importance potentially less significant. Second, just as 
is done here for region II, we have analyzed regions I and III in some detail, making use 
of the well-known forms of the Reynolds shear stress distributions in those regions. (For 
example, in region III, one has a linear distribution of the Reynolds shear stress, with some 
well- understood finite Reynolds number corrections.) It is then possible to patch together 
the results from all three regions. A summary of this work is as follows. In region I, beyond 
about y + ~ 30, dU + /dy + lies close to that implied by the classical logarithmic law. In region 
II, we have already seen that the log-law can be a good approximation (and even exact in 
a small neighborhood of Up)- In region III, dU + /dy + lies close to the classical value for 
y + smaller than 0.15.R*. It is therefore reasonable that one can smoothly match the three 
regions and obtain a logarithmic variation of velocity between y + ~ 30 and y + ~ 0.15/?*, 
regardless of whether it is exact. 

What can be said of power-laws, in particular the Barenblatt profile? In our present 
view, a self-consistent way of understanding the origin of this profile is to use Eq. (|6|) and 
the appropriate expansions for A and y p in Bl (and the appropriate version of y a in region 
B2). We believe that the rich structure of power-law profiles makes them fit the data well. 
However, a close analysis of region III ( Sahay fc Sreenivasan 1996| ) reveals that asymptoti- 



cally the velocity profile cannot be exactly logarithmic in B2. The most significant point in 
favor of power laws is that they account, in some fashion, for viscous effects to pervade in 
the classical overlap region. 

The importance of viscous effects in a region traditionally thought to be inviscid has 
other obvious analogies. The principal analogy is the Kolmogorov spectral cascade. Just 
as viscous effects are centered around the peak position of the Reynolds shear stress in 
the present problem, one may imagine the viscous effects in Kolmogorov turbulence may 
be centered around the position of maximum energy transfer in wave number space. This 
provides a natural mechanism for viscous effects to encroach the inertial region of the spectral 
space. Details are currently being worked out. 

We have so far not speculated on the physical origin of the viscous effects in the neigh- 
borhood of y p . Two possible scenarios present themselves. It may well be that the bursting 
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of wall-layer streaks (which undoubtedly extend beyond the sublayer) carry viscous effects 
with them as they penetrate up to y + = 0{Rj ). Alternatively, the viscous cores of the 
coherent vortices in the 'critical layer' region are such that their effects do not vanish at 
any Reynolds number (Barenblatt 1993). Whatever the mechanism, it appears (Sreenivasan 
1987) that the 'critical layer' or the mesolayer plays an important role in the dynamics of 
wall flows that cannot be subsumed in the classical picture. 
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